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ABSTRACT 

A searchlight type sonar is one of the systems that 
small navies use to counteract the danger which submarines 
present to their lines of supply and transport. 

In this paper, a standard search pattern for this 
type of sonar is compared with search patterns which are 
based on a consideration of the tactical value of detecting 
a submarine as a function of the relative location of the 
submarine. 

The results of the comparison suggest that is possible 
to increase the effectiveness of a searchlight type sonar 
by using a search pattern in which the sweep time allocated 
to a search sector is based on the sectors tactical importance. 
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I. 



INTRODUCTION 



Antisubmarine Warfare (ASW) training constitutes a 
major part of many navies peacetime training. One reason 
for this is that the ability to transport troops and supplies 
by sea during wartime could be a decisive factor in determining 
the outcome of the conflict. 

Several ways exist for navies to improve their ASW 
capabilities other than to obtain new ASW sensors and weapons 
systems; one of these is to develop improved tactics for use 
with their existing systems. 

An important ASW sensor system is shipboard mounted active 
sonar. Two types of active shipboard sonar systems which are 
in use today in some small navies are the searchlight type 
sonar system and the scanning type sonar system. In the 
former, a sound pulse is transmitted into and echoes are 
received from a narrow sector in a given direction. The 
area surrounding the ship must be searched in steps, sector 
by sector. In the scanning sonar, a sound pulse is transmit- 
ted and echoes are received in an omnidirectional way. The 
area surrounding the ship is searched on each pulse. 

For the searchlight type sonar, the operator has control 
of the pattern to be followed in searching the sectors 
surrounding the ship. In this paper, an investigation of 
patterns for a searchlight type sonar which tend to optimize 
detection in certain prefered sectors is made by using com- 
puter simulation techniques. 
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II. NATURE OF THE PROBLEM 



The problem of finding an optimum search pattern for a 
searchlight type detection system was addressed by Koopman 
(Ref. 1). As a measure of effectiveness, Koopman used W, 
the effective sweep width, and he concluded that when the 
instantaneous probability of detection during dt is Adt, 
assuming no previous detection, and A = A^(r)/ r 2, where 
A^(r) decreases with increasing r, the search pattern which 
maximizes the probability of detection for an infinite 
straight line encounter model consists in fixing the line of 
sight (or sonar beam) directly along the axis of abscissas, 
dividing the time equally between the right and the left 
axis. He made the following remark "It would be misleading 
to conclude that scannings should always be confined to the 
beam. In most cases, it is imperative to detect the target 
early. . . " 

In this paper, some search patterns, which consider the 
tactical value of detecting a submarine as a function of its 
relative location are investigated. 

A standard pattern for searching with a searchlight sonar 
is to start on the starboard beam and, after the first emission 
or "ping" is sent out, train the projector ten degrees toward 
the bow, send the second "ping", again train ten degrees 
toward the bow, etc. , until the projector is aimed directly 
toward the bow. From that position, the projector is then 
trained to the port beam and the above procedure is repeated. 
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The time between pings is determined by the range scale 
selected by the operator which is usually the maximum range 
for the equipment. 

In order to compare this search pattern with other 

patterns, an area around the ship was divided into sectors 

as shown in Fig. 1. The semicircular boundary of the area is 

at a distance r from the ship and it is called the maximum 

max 

range of the sonar. It will be defined more explicitly on 
page 12. 

In each sector, for the purpose of illustrating the 
approach being considered here, a somewhat arbitrary relative 
value for the detection of a submarine in the sector was assign 
ed. The value was intended to indicate the effectiveness of a 
detection in the sector in preventing a submarine from com- 
pleting its mission. The rationale supporting the relative 
values was in part as follows: a submarine detected in Sec- 

tor B is in a favorable position to be attacked, while a detec- 
tion in Sector C provides only a minimum time to attack. A 
detection in Sector A is intermediate in value to detections 
in Sector B and Sector C. 

To make the comparison between the standard pattern and 
the pattern to be considered here, no search effort was assign- 
ed aft of the beam, also a minimum range of detection for a 
sonar will be assumed. These considerations are the basis for 
defining Sector D in which no detection can occur. If an appli 
cation of the approach considered here were to be made, the 
numbers of sectors, their limits and values attached to them 
should be determined by specific tactical considerations. 
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FIGURE 1. SEARCH SECTORS WITH THEIR RELATIVE VALUES 
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A. MATHEMATICAL MODEL 



The values which have been assigned to the different 
sectors have the following order v B >v A >v c an <3 they will 
now be specified as being all positive. 

A negative V D was assigned to the event the target 
reaches Sector D, that is, the event the target enters the 
sonar search area but is not detected. 

The problem then is to try to maximize the expected value 
of the value of detection by choice of the sweep pattern. 

To obtain the detection probability on a single ping, 
a model which is outlined by Urick (Ref. 2) was used. 

The sonar equation, as given by Urick is 

(1) 10 log S/N q = SL - 2TL + TS - (NL - DI) 

where S/N q is the signal-to-noise ratio, TL is the trans- 
mission loss, TS is the target strength, SL is the source 

level, NL is the noise level and DI is the receiving directi- 
vity index. 

In simple sonar detection models, S/N q is related to 
the probability of detection p^ and the false alarm proba- 
bility p^ . If a false alarm probability is specified, 
then the probability of detection becomes a function of the 
signal-to-noise ratio alone. By using this function and 
the sonar equation, the probability of detection can be 
expressed as a function of the transmission loss, the source 
level, the target strength, the noise level and the directi- 
vity index. 
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The transmission loss is a function of the range of the 
target and for a given range r in yards it is assumed here 
to be given by 

(2) TL = 20 log r + a r 10 ^ 

where the first term represents a spherical spreading loss 
and in the second term a is in units of decibels per kiloyard 
and the second term represents absorption loss. 

Given such a relation between TL and range r of a target 
and values for SL, TS, NL and DI and a relation between the 
probability of detection p^ and the signal-to-noise ratio 
S/N q , the probability of detection p d for a single ping can 
be expressed as a function of the range to the target, that 
is the relation p d = p^(r) can found. 

To determine such a relation, the characteristics assumed 
for the searchlight type sonar were the following: an output 
of 1000 watts of acoustic power, a beam width of 16 degrees, 
a pulse length of 0.1 sec at a frequency of 10 kHz. and a 
receiver bandwidth of 500 Hz. The sonar was assumed to be 
hull mounted on a destroyer traveling at 18 knots. The 
angular width of the beam was defined as the angle between 
the 3 db. down rays. 

The sonar equation is often written as 

(3) TL = 1/2 (SL + TS - NL + DI - DT) 
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where DT = 10 log (S/N ) and (S/N ) is the signal-to- 
noise ratio required to achieve some operating point, that 
is, to achieve a particular pair of values of (Pf /P d ) • The 
significance of the equation in this form is that it relates 
the values of TL, SL, TS, NL and DI to a particular signal- 
to-noise ratio and, hence, when a relation between S/N o and 
Pf and have been specified, to a particular operating 
point. For example, if an operating point (p^/P,) is desired 
and a relation between S/N q , p f , and p^ and values of SL, 

TS , NL and DI are specified, then the value of TL which will 
give the desired operating point can be found. The value 
of p^ at the operating point is often said to be the value 
required in order that the target will be "just detectable". 

In this paper, it will be assumed that if p^ is less 
than .1, then it is in effect zero. This implies that, 
for a given p^, there is a maximum range of detection r max 
if p^ is a non-increasing function of the range, which will 
be the case for the model used here. 

A value for p^ was chosen by specifying that the proba- 
bility of more than two false alarms in one hour was .1. 

-5 

This implies the value of p^ is 1 (10) . This is shown in 

Appendix A. 

The Source Level can be expressed as (Ref. 1, p. 63) 

SL = 71.5 + 10 log P + DI t . 

For an acoustic power of 1000 w and a DI T = 20 db for the 
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central ray of the beam, SL = 118.5 db where the reference 

-12 2 

level is 0.64(10) watt/cm . From reference 1, Fig 11.10, 
the self noise level NL for a destroyer at 18 knots is -40 
db at 25 kHz. 

A noise level slope of -6 db/octave will be assumed in 
this paper. To compute the value of NL at 10 kHz, consider 
the following argument. Let x represent the number of octaves 
between 10 kHz and 25 kHz, then 

2 X = f 1 /f 2 =2.5 

and x = 1.32. Hence at 10 kHz 

NL = - 40 -(-6) (1.32) = - 32.08 db. 

The target will be assumed to be a submarine at bow 
aspect with target strength TS = + 10 db (Ref. 1, Table 9.2). 

The following model was adopted to relate signal-to-noise 
ratio to p^ and p d (Ref. 4) 

(4) p d = 1 -4> (v fc - /d ) 

(5) p f = 1 -4> (v t ) 

(6) d = wt (S/N) 2 

where w is the bandwidth, t is the pulse length and N = wN o « 

-5 

For a given value of p f which in this case is 1* (10) , v fc 

is determined by equation (5) , and with v fc and given values 
of w and t, p d became a function of S/N o through equation 
(6) and (4) . 
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This model corresponds to that for a narrow band gaussian 
signal in gaussian white noise where S/N<<1. The function 
= p^ (r) could be defined by first computing values of 
TL for various values for r by using equation (2) and then 
using equation (1) to compute values of 10 log S/M . With 
these values and equation (6) values of d could be computed 
and with the value of v^ obtained from equation (5) , the 
values of p^ which corresponded to the various values of r 
could be determined by using equation (4) . 

Because of the way the function was to be used, particular 
values of p^ were first specified and then values of r which 
corresponded to them were determined. These values, along 
with the corresponding values of TL are listed in Table 1. 

The value adopted for r is that determined by the 
operating point (1(10) .1}. 

Since the sonar equipment of the type being considered 
here is usually limited in its detection capability for 
targets at very close range, a minimum detection range r 
of 200 yards was adopted in the model. 

To simplify the problem, it was assumed that the target 

was a point target and that if a target is in the beam for 

m pings, the probability of detecting it is 

m 

P {Detection} = 1 -^J 1 ^ { 1 - 

where both m and r^ are determined by the relative track 
of the target through the area scanned by the sonar beam. 

As a further simplification, the targets track was assumed 
to be parallel to the ships track, so that a target will enter 
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table 1 




Table 1. Values of the transmission loss and ran9e 
for particular values of the probability of detection 



and a false alarm probability 
of more than two false alarms 
equal to . 1 . 



such that the probability 
in a one hour period is 
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the semicircular area swept by the beam if its lateral range 

x is such that |x|<r 

1 1 max 

In order to determine the number of scans m on a target 
during its straight line encounter with the searcher, the 
sonar position in its search cycle at the moment the target 
crosses the semicircular boundary must be specified. In 
addition, the targets lateral range, sonar scanning period 
and the sweep pattern must be specified. It is reasonable 
to assume that given the target enters the search area, its 
lateral range is a random variable which is uniformly dis- 
tributed between -r and r and that the target crosses 

max max 3 

the line which is tangent to the semicircular area swept by 
the beam and parallel with the x-axis at a time which is 
uniformally distributed over the time for a complete scan 
cycle. 

Even with the above assumptions, the problem of determin- 
ing the average of the value for a pattern is difficult to 
handle analytically. However, a Monte Carlo simulation can 
provide a satisfactory solution to this problem. 

In this investigation, a computer was used to do a 
Monte Carlo simulation in order to estimate the expected 
value of the value for some patterns. 

This simulation and its results are discussed in the 
remainder of the paper. 
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III. DESCRIPTION OF THE SIMULATION 



A search pattern which might be a realizable optimum 
one under the assumed conditions can be described as follows: 
start with one normal sweep beam to beam, and then allocate 
one partial sweep between the rays which bound Sector B. 

This will be called a 2/1 search pattern since Sector B is 
searched twice for each search of Sector A. Using this 
terminology, the standard search pattern would be called a 
1/1 search pattern. By allocating additional partial sweeps 
to Sector B, the search patterns 3/1, 4/1, etc., are generated. 

It has been assumed that a submarine can enter the sonar 
area of detection at any position of the projector and at 
any lateral range. To simulate this, a target was generated 
at random on a line parallel to the x-axis and tangent to 
the semicircular area swept by the beam each time the projector 
was advanced to a new position. 

To reduce starting transient effects and to allow each 
target an opportunity to be detected, the search process 
was not started until the first target generated had advanced 
approximately 1750 yards and the last target generated was 
allowed to travel the same distance before the search process 
stopped. 

The simulation was run for the search patterns 1/1, 2/1. 
3/1, 4/1, and 5/1 with 500 targets generated for each pattern. 
Ten different runs of each pattern were made in order to 
obtain a statistically adequate sample. 
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The sectors chosen and the detection values assigned to 
them were as follows, where bearings are relative to the 
heading of the ship: Sector A: between bearings 040° and 

090° and between 270° and 320° and between ranges of 800 
and 1800 yards. Value: 2. Sector B: between bearings 

320° and 040° and between ranges of 800 and 1800 yards. 

Value: 10. Sector C: between bearings 270° and 090° and 

between ranges of 200 and 800 yards. Value: 1. Sector D: 

between bearings 270° and 090° and between ranges of 0 to 
200 yards and all the area aft of the x-axis. Value: -1. 

The geometry used for the simulation is shown in Fig. 2. 
The ship is at the origin of a system of rectangular coordi- 
nates which is moving in the positive direction of the y-axis 
at a speed of 18 knots. 

As only positive values of x need to be considered because 
of the symmetry of the sectors and sweep patterns, the area 
of interest for the simulated search is the positive quad- 
rant and targets were uniformly generated on the positive 
section of the line parallel to the x-axis. This section 
of the line is 1800 yards long and it is 1800 yards from 
the ship. The value 1800 was used since it is a little 
larger than the range corresponding to p d = 0.1 which was 
arbitrarily said to be equivalent to p d = 0. 

The targets' moved in the negative direction of the y- 
axis and parallel to it at a speed of 10 knots, so that their 
relative speed during an encounter was 28 knots. Both ships 
kept their courses and speed constant. 
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FIGURE 2. GEOMETRY OF THE SIMULATION 
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The name of the variables and data used are given in 
Appendix B and the schematic flow chart is shown in Fig. 3. 

The matrix TAR and vector ITAR were used to keep the 
parameters of each target generated. The coordinates of 
the position of the target were kept in TAR. The condition 
of the target was said to be active if it was in the detec- 
tion area and not detected and passive if it was in the 
detection area and detected or if it had crossed into 
Sector D undetected. 

When the search process starts, the first ping is on 
the starboard beam, for each ping after the first ping, the 
projector is trained 10 degrees and the target advanced 
downward 35 yards which corresponds to the relative distance 
traveled by the target in a time interval equal to that 
necessary for a sound pulse to go to 1800 yards and return. 

A sound speed of 4800 ft/sec was used. 

On each ping, the angular positions of the active targets 
were computed and compared in order to determine those, if 
any, which were inside the sonar beam. A beam width of 
16 degrees was used. All of those targets which were out- 
side of the beam could not be detected and they were advanced 
downward. If a target was inside the beam, its range to the 
ship was computed. If this range were less than the minimum 
range of 200 yards, no detection was possible and the target 
was then assigned to Sector D. 

For targets whose ranges were greater than the minimum, 
the transmission loss to the target was computed and the 
probability of detection was obtained. If this probability 
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was less than 0.1, no detection occurred and the target 
was advanced. 

To determine if a detection event occurred, a random 
number, uniformly distributed (0,1), was generated and if 
the random number was less than the probability of detection 
a detection occurred. Otherwise a detection did not occur. 
Detections were assigned to the sectors in which they 
occurred . 

After the above process was completed, the projector 
was trained, a new ping emitted and the process was started 
again. 

Whenever the projector reached the position in which it 
was trained toward the bow of the ship, the search on the 
port side was simulated. This was done by executing the 
generation and advance of the targets in the positive quad- 
rant for the numbers of pings required to cover the port 
side for the pattern being used. Detection of targets could 
not be made during the period necessary to generate these 
pings. The above procedure was possible because of the port/ 
starboard symmetry of the problem. 

The targets which crossed the area without being detected 
were assigned to Sector D. 

The search pattern execution was controlled by the part 
of the program labelled Control of the Search Sweep Pattern. 
It essentially counted the times the projector went into 
each sector and produced the desired sweep ratio. This 
section also controlled the simulated search on the port 
side. 
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When the five search patterns had been realized, a new 
run for all the sweep patterns was executed. 

The detection values obtained for each sweep pattern 
are kept in the matrix STORE. This matrix and the average 
of the value and standard deviation of the value for each 
pattern was the output of the simulation. 

The outputs of a given search pattern for the ten 
different runs can be considered to be samples from the 
population of values for that sweep pattern. 

The Mann-Whitney test was performed with the data in 
the matrix STORE and the result was printed. 
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FIGURE 3. SCHEMATIC FLOW CHART (5) 
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IV. CONCLUSIONS 



The averages of the values of detection for the five 
sweep patterns investigated and for the assigned values 
suggested in some cases that the expected value of the value 
of detection is different than that of the 1/1 standard 
pattern. To test this conjecture, the Mann-Whitney test 
was used. This test was used because it is nonparametric , 
the distribution of the value of detection was not known, 
and the sample size was not large enough to assume a normal 
distribution . 

The Mann-Whitney test (Ref. 3) with the assumption that 
if there is a difference between population distribution 
functions, then that difference is a difference in the means 
of the distributions allows one to make the following hypo- 
thesis: 

H q : E{X (1/1 ) } 1 E{ Y (n/1) } 

H 1 : E {X (1/1 ) } < E { Y (n/1) } 

where X(l/1) is the random variable representing a value 
from the population of values of detection with the (1/1) 
standard search pattern and Y(n/1) is the random variable 
representing a value from the population of values of detec- 
tion obtained with a (n/1) search pattern. 

If the null hypothesis is accepted, that is if the 
expected value of X is greater than or equal to the expected 
value of Y, this means that at the level of significance 
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chosen, there is not a significant improvement in the value 
of detection with that pattern. If the null hypothesis is 
rejected, this means that there is a significant improve- 
ment by using that sweep pattern. 

For the test, a 90% confidence level was chosen. 

The Mann-Whitney test results are given in Table II 
below. 

TABLE II 

RESULT OF MANN-WHITNEY TEST 



VALUE OF T STATISTIC 50.00 

FOR SWEEP PATTERN 1/1 ACCEPT H Q 

VALUE OF T STATISTIC 43.50 

FOR SWEEP PATTERN 2/1 ACCEPT H Q 

VALUE OF T STATISTIC 36.00 

FOR SWEEP PATTERN 3/1 ACCEPT H Q 

VALUE OF T STATISTIC 32.00 

FOR SWEEP PATTERN 4/1 REJECT H Q 

VALUE OF T STATISTIC 34.50 

FOR SWEEP PATTERN 5/1 ACCEPT H Q 



It should be noted again that the above results were obtained 
with arbitrary sectors and arbitrary sectors values of 
detection and that the encounters were particular straight 
line encounters. In addition, the assumed characteristics 
of the sonar equipment was specific and the relation between 
TL and r was restrictive. 
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Estimates of the probability that a target which enters 
a search pattern will be detected can be made as follows: 
Assign the values = Vg = = 1 and V D = 0 and then 

divide the average of the values which are obtained from 
the simulation by the number of targets generated in the 
simulation. The resulting numbers are the estimates. 
Estimates for the five sweep patterns investigated are 
given in Table III. 



TABLE III 



ESTIMATES OF THE PROBABILITY OF DETECTION 



Sv/eep Pattern 



Est. Prob. of Detection 



1/1 



0.91 



2/1 



0.89 



3/1 



0.85 



4/1 



0.82 



5/1 



0.79 
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APPENDIX A 



Determination of False Alarm Probability. 

Since p^ and p^ apply to a resolution cell and the 
resolution cell's size in time is equal to the pulse length, 
in one hour n = 3600 / .1 = 36000 resolutions cells will be 
sampled. 

If independence is assumed and if X represents the 
number of false alarm in one hour period then the probabi- 
lity of less than- two false alarms in one hour is 



P [X<2] 





n-i + 




Since p^<<l and n>>l 



P [X<2] » (1-P f ) n + n p f (l-p f ) n 



(l-p f ) n (1+n p f } 



(l-np f ) (l-np f ) 



2 

Hence, p^ 



1 - P [ X <.2 ] 

2 



and if P[X<2] = .9, that is 



n 



P [X>2] - .1 then 



p f “ 0.8784.10 



-5 



or 



p f “ 1.10 



-5 
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APPENDIX B 



DATA 


VARIABLES AND DATA IN THE PROGRAM 


NSS 


= Number of targets generated 


RMAX 


= Maximum sonar search range 


ATT 


= Attenuation factor for propagation 


LUPA 


= Upper limit of Sector A (starboard side) 


LLUPA 


= Upper limit of Sector A (port side) 


IX 


= Seed for random number generator 


LIM 


= Numbers of additional targets generated 



(LIM=NSS+50) 



YS 


= y distance of generation of target 


RMIN 


= Minimum sonar detection range 


RDA 


= Upper range limit of Sector C 


VSH 


= Speed of the ship 


VSS 


= Speed of the submarine (target) 


VSO 


= Speed of sound 


ITETA2 


= Half beam width 


INCR 


= Step in train of the projector 


NUM 


= Number of runs 


JN 


= Increment in number of ping. 


TTAB 


= Value of table for Mann-Whitney test 


WEn 


= Value of detection in Sector n (A,B,C,D) 


MOSn 


= Required number of sweep in Sector n (A,B) 


TLxx 


= Transmission loss corresponding to O.xx prob 




of detection 
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VARIABLES 



I BETA 


= Angular position of projector 


IBEUB 


= Upper ray of sonar beam 


IBELO 


= Lower ray of sonar beam 


ALFA 


= Angular position of the target 


ADV 


= Advance of target in each inter ping time 


PD 


= Probability of detection 


YFL 


= Random number U (0.1) 


TL 


= Actual transmission loss 


RT 


= Range of target 


M 


= Number of targets actually generated 


IGEN 


= Counter for number of pings 


T 


= Time between pings 


TVTEn 


= Total value of detection in Sector n 


TTW 


= Total value of detection 


NOn 


= Number of detections in Sector n 


NOSn 


= Counter to control sweep pattern 


T 


= Statistic of Mann-Whitney test 



MATRICES AND VECTORS 



STORE (I,J) 


= Value of detection for run I and sweep 
pattern J 


VALUE (J) 


= Average value of detection for sweep 
pattern J 


SIGMA (J) 


= Standard deviation for value of detection 
of pattern J 


TAR (K,L) 


= Data of each target K 
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* 









ITAR (K) 



A (I ) 
R(I) 



(K , 1 ) x coordinate of target K 
(K,2) y coordinate of target K 
= Condition for each target K 

f 

1 = Target in the area undetected (Active) 
< 0 = Target detected or crossed the area 
undetected (Passive) 

= Data of run (I) for Mann-Whitney test 
= Working vector for Mann-Whitney test 



SUBROUTINES 

RANK =. Subroutine IBM to rank a vector of values 

RANDU = Subroutine IBM to generate a random number 

U(0.1) 



32 



COMPUTER OUTPUT 



MATRIX OF VALUES OF DETECTION 



RUNS 




SWEEP 


PATTERNS 








1/1 


2/1 


3/1 


4/1 


5/ 1 


1 


2959.00 


2590.00 


2684. 00 


2825.00 


2812.00 


2 


2748. 00 


2760. 00 


2831. 00 


2837.00 


2766.00 


3 


2853.00 


2786.00 


2823.00 


2713.00 


2781.00 


4 


2780. 00 


2799. 00 


2791.00 


2678.00 


2811.00 


5 


2587.00 


2707 .00 


2934.00 


2849.00 


2866.00 


6 


2590.00 


2684.00 


2825. 00 


2812.00 


2775. 00 


7 


2760. 00 


2831.00 


2837. 00 


2776.00 


2788.00 


8 


2786.00 


2823.00 


2713.00 


2871.00 


2806. 00 


9 


2779.00 


2791.00 


2678.00 


2811.00 


2566. 00 


10 


2707. OO 


2934.00 


2849.00 


2886.00 


2779. 00 


S 


. PATT 


AV. VALUE 


STD. DEV 








1/1 


2754.899 


111.160 








2/1 


2768.499 


93. 123 








3/1 


2796.499 


81.451 








4/1 


2805. 799 


66.491 








5/1 


2786.999 


86.660 








RESULT OF 


MANN-WHITNEY TEST 







VALUE OF T STATISTIC 50.00 
FOR S. PATTERN 1/1 ACCEPT HO 

VALUE OF T STATISTIC 43.50 
FOR S. PATTERN 2/1 ACCEPT HO 

VALUE OF T STATISTIC 36.00 
FOR S. PATTERN 3/1 ACCEPT HO 

VALUE OF T STATISTIC 32.00 
FOR S. PATTERN 4/1 REJECT HO 

VALUE OF T STATISTIC 34.50 
FOR S. PATTERN 5/1 ACCEPT HO 
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COMPUTER PROGRAM 



C SIMULATION GF SEARCH AND DETECTION OF A 

C SEARCHLIGTH TYPE SONAR 



DIMENSION TAR{ 550, 2), ITARI550) 

DIMENSION STORE (5,10), 1 0( 5 ) , V ALUE ( 5 ) , S I GM A { 5 ) 
DIMENSION A( 23 ) , R( 20) 

DATA TAR/ 1100*0 .0/ , I TAR/ 5 5 0*0/, I 0/ 1 , 2, 3 , 4 , 5/ 

DATA ST ORE/ 5 0* 3. 0/ , VALUE/5*0.0/, SIGMA/ 5*0.0/ 

DATA INCR/10/ , J N / 1 / , ITETA2/8/, ATT/7.0/ ,MOSA/l/ 

DATA RMAX/1800.0/ , RD A/ 800. 0/,RM IN/200.0/ 

DATA VSH/18.0/,VSS/10. 0/,V SO /16 00.0/, YS/18 00.0/ 
DATA TL10/77.1 1/,TL 20/ 76.81 /,TL 30/76. 62/ 

DATA TL40/76. 4 7/ , T L 5 0/ 76 . 3 4/ , T L6 0/76.21/ 

DATA TL70/76.09/, TL80/75.95/ ,TL90/75.77/ 

DATA WE A/2 .0/ , WEB/ 1 0.0/ , W EC/ 1. 0/ , WED/- 1 . 0/ 

NSS=500 
NUM= 1 0 
LUP A = 5 0 

LLUPA=180-LUPA 
T=RMAX*2/VS0 
LI M=NSS+50 
LUM=NUM- 1 

C GENERATE RUNS 

DO 7000 KK = 1,\IUM 
READI5 ,8100 ) IX 

C GENERATE ALL TARGETS FOR EACH SWEEP PATTERN 

DO 2000 J J =1 ,5 

M0S8 = J J 

I BET A= 0 

NO SA = 0 

N0SB=0 

NOA=0 

NOB = 0 

N0C = 0 

NOD=0 

M=0 

I GE N = 1 

C TARGETS ARE GENERATED AT ALL POSITIONS OF THE BEAM 

180 M=M+ 1 

IF(M.GT.LIM) 30 TO 1900 
IF(M.GT.NSS) 30 TO 210 

C GENERATE TARGET UNIFORMLY DISTRIBUTED 

CALL RANDU (IX, I Y , YF L ) 

IX= I Y 

TARIM, 1 ) =YS* YF L 
T AR ( M, 2 ) = YS 
I T AR ( M ) = 1 

THE SEARCH START FIRST TARGET HAS ADVANCED 50 STEPS 
APROX 1750 YARDS 

IF ( I GEN- 50} 1,210,210 

1 I BET A = 0 
GO TO 150 
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C CONTROL OF THE SEARCH SWEEP PATTERN 

210 I F ( ( NOSA.EQ.MOSA) .AND. (NOSB. EQ. MOS B) ) GO TO 265 
I F ( IBETA.LE.9D ) GO TO 200 

IF ( ( IBETA. GT .90 ) .AND. ( I BETA.LE .LLUPA ) ) GO TO 150 
I F ( IBETA.GT.180; GO TO 260 
IFJNOSA.LT. MOSA) GO TO 150 

C SOME SECTOR NEED EXTRA SWEEP 

I BETA = LUPA 
NO SB = N0 SB+ 1 

IF(NOSB.EO.MOSB) I B ET A= 0 
GO TO 200 

C BOTH SWEEP COMPLETED, START A NEW CYCLE 

265 NO SA = 0 
NOSB = 0 
GO TO 200 

C BOTH SECTOR NEED EXTRA SWEEP 

260 NOSA=NOS A+ 1 
N0SB=N0S3+ 1 
I B ET A = 0 

IF(NOSA.EQ.MOSA) I BETA=LUPA 

IFJ (NOSA.EQ.MOSA). AND. ( NO SB . EQ . MOSB ) ) I B ET A = 0 

C DETERMINE POSITION OF THE SONAR BEAM 

200 IBEJP=IBETA+ITETA2 
I B E L0 = I BET A-I T ETA2 
IFJ I3EL0.LT. 0) I BE L C =0 

CHECK IF TARGET IS IN THE BEAM 
CHECK ALL ACTIVE TARGETS 

DO 120 J = 1 » M 

TARGET CROSSED AREA UNDETECTED 
ASSIGN TO SECTOR D 

IFJITARJJ)) 120,120,101 
101 IFJ TARJ J ,2 ) .GT .30.0) GO TO 110 
GO TO 258 

110 ALFA=AT ANJ TARJ J,2)/TAR(J,1) ) *360 . 0/6. 2832 
220 IFJ (ALFA. LE. I3EUP) .AND. (ALFA. GE. IBELO) ) GO TO 300 

C ADVANCE TARGET 

250 AD V = T* { VSH + VS S) *2 OOO. 0/3600.0 
TARJ J, 2 )=TAR( J ,2)-ADV 
IF (TARJ J,2 ) .GT .30.0) GO TO 120 
GO TO 258 

C IF TARGET IS IN THE BEAM COMPUTE RANGE 
3 00 RT = SORT ( TAR { J , 1 ) ** 2+T AR ( J , 2 ) **2 ) 

C NO DETECTION FOR MINIM. RANGE . ASSIG TO SECTOR D 
IFJRT.LT. RMIN)GO TO 258 
C COMPJTE TRANSM I SION LOSS 

TL=20.0*AL0G1D ( RT ) +ATT *RT/ 1 000 . 0 
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C CHECK TL TO FIND PROB. OF DETEC. 

C ASUME NO DETECTION IF PROB. DETEC IS 0.1 OR LESS 



410 
42 0 
430 
440 
450 
46 0 
480 
49 0 



47 0 



I F ( TL. GT.TL 10) GO TO 250 

IFITL.LT. TL90) GO TO 410 

I F ( TL.LT.TL30) GO TO 420 

IFITL.LT. TL70) GO TO 430 

IFITL.LT. TL60) GO TO 440 

IF! TL . LT. TL 50) GO TO 450 

IFITL.LT. TL40) GO TO 460 

I F ( TL. LT.TL30) GO TO 480 

IFITL.LT. TL20) GO TO 490 

PD=0 . 1 

GO TO 470 

P D=0 . 9 

GO TO 470 

PD=0. 8 

GO TO 470 

PD=0 . 7 

GO TO 470 

P D=0 . 6 

GO TO 470 

PD=0. 5 

GO TO 470 

PD =0 . 4 

GO TO 470 

PD=0 .3 

GO TO 470 

PD= 0. 2 

GENERATE RANDOM NUMBER UI0,1) TO DECIDE 

CALL RANDU I IX, I Y , YF L ) 

IX = I Y 

IFIYFL.GT. PD) GO TO 250 

EVENT DETECTION - ASSIGN TO ONE SECTOR 



600 IFIRT.LT. RDA) GO 
I F{ ALFA .LE .LUPA ) 
; SECTOR B 

N0B=NG3 + 1 
GO TO 640 
; SECTOR C 
615 NQC = NOC + 1 
GO TO 640 
; SECTOR A 

620 NOA=NOA+ 1 
GO TO 640 



TO 

GO 



615 

TO 620 



SECTOR D 
NOD=NOD+1 
ITARI J) =0 
CONT INUE 
GO TO 160 

WHEN BEAM IS AT PORT SIDE, ADVANCE TARGETS 

DO 130 J=1 , M 
IFIITARIJ)) 130,130,140 
ADV=T*( VSH+VSS ) -20 00 . 0/3600 . 0 
TARI J,2)=TAR( J.2J-ADV 
TARGET CROSSED AREA NO DETECTION 
IF ( TAR ( J,2).GT.30.0) GO TO 130 
NOD=NOO+ 1 
IT AR I J ) =0 
130 CONT I N J E 



258 

640 

120 



150 

140 



DETECTION 



I N STBD. SIDE 
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C ROTATE BEAM 

160 IGEN=IGEN+JN 

IBETA=I BETA+IMCR 
GO TO 180 

C COMPUTE VALUE OF DETECTION AND STORE IT 

1900 TWEA=WEA*NOA 
T WEB=WEB~NOB 
TWEC=W E C-NOC 
TWED=WED*NOD 
TTW=TWEA+TWEBf TWEC+TWED 
STOREKK, JJ)=TTW 
2000 CONTINUE 
7000 CONTINUE 

C COMPUTE MEAN AND STANDARD DEVIATION 

DO 6000 J = l,5 
DO 5000 K=1,NJM 

VALUE! J ) =VALUE ( J )+ STORE! K, J) /NUM 
5000 CONT IN J E 

DO 6550 I K = 1 , NUM 

SIGMA! J ) =S I GM A ( J ) + ( STORE ( I K, J) -VALUE ! J) ) **2 
6550 CONT IN J E 

SIGMA! J ) =SQRT! S I GMA { J ) / LUM ) 

6000 CONTINUE 

C PRINT MATRIX OF VALUE OF DETECTION 

WR IT E ( 6 , 9999 ) 

WRITE l 6,9990) 

WRITE! 6, 8500) ( 10! J ) ,U = 1 ,5) 

DO 8000 LL= 1 , NJM 

WRITE (6,8200) L L , l STO RE ( L L , JK ) , JK=1 , 5 ) 

8000 CONTINUE 

C PRINT MEAN AND STANDARD DEVIATION OF DETECTION VALUES 

WRITE (6,8300) 

DO 6500 U = 1 , 5 

WRITE (6,8400) I 0 ( U ) , VAL UE ( J ) , SI GMA { J ) 

6500 CONTINUE 

C MANN - WHITNEY TEST 

HO: E ( VALUE 1/1) GREATER OR EQUAL TO E ( VALUE X/l) 
Hi: E ( VA LUE 1/1) LESS THAN ElVALUE X/l) 

C FOR 90 % OF CONFIDENCE (ONE TAIL TEST) 

N 1 = 1 0 

N2=10 

N = N 1 + N 2 

TT AB=3 3 

WRI TE ( 6 ,9000 ) 

C READ VECTOR SWEEP PATTERN 1/1 

DO 15 J = 1 , N1 
A( J )=STORE ( J, 1 ) 

15 CONTINUE 
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C READ VECTOR TO TEST (FIRST IS PATTERN I/I) 

N2=N2+ 1 
DO 20 1=1,5 
DO 3 3 I J =N2 , N 
JI=IJ-10 

A(IJ)=STORE(JI , I ) 

33 CONTINUE 

C RANK OF VALUES 

CALL RANK( A,R, Nl ) 

C SUM RANK OF FIRST VECTOR 

R2=0. 0 
DO 10 M=1,N1 
R2=R2+R ( M) 

10 CONTINUE 

C COMPUTE T STATISTIC 

T = R2-N1*((N1 + 1 ) / 2 . 0 ) 

WR I TE ( 6 » 93 1 3 ) T 
IF(T.GE.TTAb) SO TO 50 
WRITEtS, 9100)1 
GO TO 23 

50 WRITE{ 6, 9200) I 
20 CONTINUE ' 

C FORMATS 

8103 FORMAT (19) 

8200 FORMAT ( IX , // , 14X, 14, 5F9. 2 ) 

8300 FORMAT (' 1' ,//// ,20X, 'S. PATTERN AV. VALUE S.DEV') 
8433 F0RMAT(1X,//,24X,I2,3X,2F13.3) 

8500 FORMAT ( IX, // , 16X, 519) 

9000 FORMAT ( ' 1' ,/////, 25X,' RESULT OF MANN-WHITNEY TEST') 
9313 FORMAT ( IX ,/////, 2 1 X ,' VALUE OF T S T AT I ST I C ' , FI 0 . 2 ) 
9100 FORMAT ( IX, /, 21 X, ' C 0R S . PATTERN • , I 2 . • /I REJECT HO') 
9200 FORMAT ( IX, / ,21 X, ’ FOR S . PATTERN ', 12 ,'/ 1 ACCEPT H3 ' ) 
9990 FORMAT ( IX,///, 14X, 'RUNS' ,15X,' SWEEP PATTERN') 

9999 FORMAT (' 1' ,/////, 28X, 'MATRIX OF VALUE OF DETECTION') 
STOP 
END 



} 
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